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INTRODUCTION 

This application note describes the implementation of 
various digital filters using the PIC17C42, the first 
member of Microchip's 2nd generation of 8-bit 
microcontrollers. The PIC17C42 is a very high speed 
8-bit microcontroller with an instruction cycle time of 
250 ns (@ 16 MHz input clock). Even though the 
PIC17C42 is an 8-bit device, its high speed and effi- 
cient instruction set allows implementation of digital fil- 
ters for practical applications. Traditionally digital filters 
have been implemented using expensive Digital Signal 
Processors (DSPs). In a system the DSP is normally a 
slave processor being controlled by either an 8-bit or 
16-bit microcontroller. Where sampling rates are not 
high (especially in mechanical control systems), a sin- 
gle chip solution is possible using the PIC17C42. 

This application note provides a few examples of imple- 
menting digital filters. Example code for 2nd order Infi- 
nite Impulse Response (IIR) filters is given. The 
following type of filters are implemented: 

• Low Pass 

• High Pass 

• Band Pass 

• Band Stop (notch) filter 

This application note does not explain how to design a 
filter. Filter design theory is well established and is beyond 
the scope of this application note. It is assumed that a filter 
is designed according to the desired specifications. The 
desired digital filters may be designed using either stan- 
dard techniques or using commonly available digital filter 
design software packages. 

Finite Impulse Response (FIR) filters have many 
advantages over IIR filters, but are much more resource 
intensive (both in terms of execution time and RAM). On 
the other hand, IIR filters are quite attractive for 
implementing with the PIC17C42 resources. Especially 
where phase information is not so important, IIR filters 
are a good choice (FIR filters have a linear phase 
response). Of the various forms used for realizing digital 
filters (like, Direct form, Direct II form, Cascade form, 
Parallel, Lattice structure, etc.) the Direct II form is used 
in this application note. It is easy to understand and 
simple macros can be built using these structures. 



THEORY OF OPERATION 

Digital filters in most cases assume the following form 
of relationship between the output and input 
sequences. 



M N 
y(n) = £ a.y(n-i)+ £ bjx(n-j) 
i = o j = o 



The above equation basically states that the present 
output is a weighted sum of the past inputs and past 
outputs. In case of FIR filters, the weighted constants 
ai = and in case of IIR filters, at least one of the ai 
constants is non zero. In case of IIR, the above formula 
may be rewritten in terms of Z transform as: 



H(z) ~m~ 



M 
k = o 



N 

1+ X a k Z 
k = l * 



,-k 



The above equation can further be rewritten in 
difference equation format as follows: 



M N 
y(n) = £ a.y(n-i)+ £ bjx(n-j) 

i= 1 j=o 
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Realization of the above equation is called the Direct 
Form II structure. For example, in case of a second 
order structure, M=N=2, gives the following differ- 
ence equations: 

EQUATION 1: 

d{n)= x(n) + ajd(n - 1) + a 2 d(n - 2) 

EQUATION 2: 

y{n)= b d(n) + bjd(n-l)+b 2 (d(n-2)) 

The above difference equations may be represented as 
shown in Figure 1 . 



FIGURE 1 : 2ND ORDER DIRECT FORM II 
STRUCTURE (TRANSPOSED) 
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The structure as shown in Figure 1 may be cascaded to 
attain a higher order filter. For example, if two stages 
are cascaded together, a 4th Order MR Filter is 
obtained. This way, the output of the 1st stage becomes 
the input to the second stage. Multiple order filters are 
thus implemented by cascading a 2nd order filter struc- 
ture as shown in Figure 1 . 

IMPLEMENTATION 

A 4th order MR Filter is implemented by cascading two 
of the structures shown in Figure 1 . The output Y (out- 
put of each filter stage) is computed by direct imple- 
mentation of Equation 1 and Equation 2. Since each 
stage is similar algorithmically, it is implemented as a 
macro using Microchip's, Assembler/Linker for 
PIC17C42. This Macro (labelled B I quad) is called 
twice for implementing a 4th order filter. The output of 
the 1st stage is directly fed to the input of the second 
stage without any scaling. 

Scaling may be required depending on the particular 
application. The user can modify the code very easily 
without any penalty on speed. Also, saturation 
arithmetic is not used. Overflows can be avoided by lim- 
iting the input sequence amplitude. All numbers are 
assumed to be 1 6 bits in Q1 5 format (1 5 decimal points, 
MSb is sign bit). Thus the user must scale and sign 
extend the input sequence accordingly. For example, if 
the input is from a 12-bit AID converter, the user must 
sign extend the 12-bit input if bit 1 1 is a one. 

The B I quad macro is a generic macro and can be 
used for all MR filters whether it is Low Pass, High Pass, 
Band Pass or Band Stop. A general purpose 16x16 



multiplier routine is also provided. This routine is 
implemented as a straight line code for 
speed considerations. 

The 4th order MR filter implemented is a Low Pass Filter 
with the specifications shown in Table 1 . 



TABLE 1 : FILTER CONSTANTS 





BAND1 


BAND2 


Lower Band Edge 


0.0 


600 Hz 


Upper Band Edge 


500 Hz 


1 kHz 


Nominal Gain 


1.0 


0.0 


Nominal Ripple 


0.01 


0.05 


Maximum Ripple 


0.00906 


0.04601 


Ripple in dB 


0.07830 


-26.75 


Sampling Frequency = 2 kHz 



The Low Pass Filter is designed using a digital filter 
design package (DFDP™ by Atlanta Signal Processors 
Inc.). The filter package produces filter constants of the 
structure shown in Table 1. Table 2 shows the filter 
co-efficients that are obtained for the above Low Pass 
filter specification. 



TABLE 2: FILTER CO-EFFICIENTS 





Co-efficients 


CO 


a1 


a2 


b0 


b1 


b2 


1 


-0.133331 


0.167145 


0.285431 


0.462921 


0.285431 


2 


0.147827 


0.765900 


0.698273 


0.499908 


0.698273 



The above filter co-efficients (5 per stage) are quan- 
tized to Q15 format (i.e they are multiplied by 32768) 
and saved in program memory (starting at label 
coef f_lpass). The constants for both the stages 
are read into data memory using tlrd and tablrd 
instructions in the Initialization Routine (labelled 
initFilter). The user may read the coefficients of 
only one stage at a time and save some RAM at the 
expense of speed. 

The sample 4th order Low Pass IIR Filter is tested by 
analyzing the impulse response of the filter. An impulse 
signal is fed as input to the filter. This is simulated by 
forcing the input to the filter by a large quantity (say 
7F00h) on the first input sample, and the all zeros from 
the 2nd sample onwards. The output sequence is the 
filter's impulse response and is captured into the 
PICMASTER's (Microchip's Universal In-Circuit 
Emulator) real-time trace buffer. This captured 
data from PICMASTER is saved to file and analyzed. 
Analysis was done using MathCad™ for Windows® and 
a DSP Analysis program from Burr-Brown 
(DSPLAY™). The Fourier Transform of this impulse 
response of the filter should display the filter's fre- 
quency response, in this case being a Low Pass type. 
The plots of the impulse response and the frequency 
response are shown in Figure 2, Figure 3 and Figure 4. 

MathCad is a trademark of MathSoft, Inc. 

DSPLAY is a trademark of Burr-Brown 

DFDP is a trademark of Atlanta Signal Processing Inc. 

Windows is a registered trademark of Microsoft Corporation. 
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FIGURE 2: 
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FIGURE 3: SPECTRUM COMPUTED 

FROM IMPULSE RESPONSE 
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PERFORMANCE 

The resource requirements for filter implementations 
using a PIC17C42 is given in Table 3. These numbers 
can be used to determine whether a higher order filter 
can be executed in real-time. The same information 
may be used to determine the highest sampling rate 
possible. 

FILTER APPLICATIONS 

Digital filters find applications in many areas especially 
those involving processing of real world signals. In 
some applications like ABS systems in an automobile, 
digital filtering becomes a must. In this case elimination 
of noise (especially glitches and false readings of 
sensors) is very critical and thus becomes a require- 
ment of digital signal processing. 



TABLE 3: RESOURCE REQUIREMENTS 



Timing (Cycles) ' 


#of Filter Stages*775 + 16 


Program Memory (locations) ' 


#of Filter Stages*68 + 290 


RAM (File Registers) 


#of Filter Stages*16 + 16 



Note 1 : The above numbers do not include the 
initialization routine. 



FIGURE 4: LOG MAGNITUDE 

SPECTRUM OF IMPULSE 
RESPONSE 



Frequency Response 




64 128 192 256 320 364 448 512 
Frequency /51 2 (kHz) 



Digital filters are also needed in Process Control where 
notch filters and low pass filters are desired because 
the signals from sensors are transmitted over long 
lines, especially in a very noisy environment. In these 
cases, typically a notch filter (centering 50 Hz or 60 Hz) 
is used. In cases of eliminating background noise, a 
band stop filter (e.g., 40 Hz to 120 Hz) is used. The 
sample code given in this application note can be used 
to design a feedback control system's digital 
compensator. For example, atypical DC Motor's digital 
compensator (like a dead-beat compensator) is of 
second order and has the same filter structure that is 
implemented in this application note. 
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TABLE 4: RESOURCE REQUIREMENTS 



Filter Order 


Cycles 


Real Time 
(@16 MHz) 


Maximum 
Sampling 


Program 
Memory* 1 ) 


RAM 


2 


791 


197.75 ns 


5.05 kHz 


358 


32 


4 


1566 


391.5 us 


2.55 kHz 


426 


48 


6 


2341 


585.25 ns 


1.7 kHz 


494 


64 


8 


3116 


779.0 ns 


1.28 kHz 


562 


80 


10 


3891 


972.75 ns 


1.0 kHz 


630 


94 



Note 1 : The above numbers do not include the initialization routine. 
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Please check the Microchip BBS for the latest version of the source code. Microchip's Worldwide Web Address: 
www.microchip.com; Bulletin Board Support: MCHIPBBS using CompuServe® (CompuServe membership not 
required). 



APPENDIX A: iir.lsOt 

MPASM 01.40 Released IIR.ASM 1-16-1997 14:48:37 PAGE 1 



LOC OBJECT CODE 
VALUE 



LINE SOURCE TEXT 



00000001 
00000000 
00000000 
00000001 



00001 
00002 
00003 
00004 
00005 
00001 
00002 
002 64 
00006 
00007 
00008 
00009 
00010 
00011 
00012 
00013 
00014 
00015 
00016 
00017 
00018 
00019 
00020 
00945 
00021 
00022 
00001 
00002 
00003 
00004 
00005 
00006 
00007 
00008 
00009 
00010 
00011 
00012 
00013 
00014 
00015 
00016 
00017 
00018 
00019 
00020 
00021 
00022 
00023 
00024 
00025 
00026 
00027 
00028 



TITLE "Digital IIR Filter Using PIC17C42" 
LIST P=17C42, columns=120,WRAP, R = DEC 
#include <pl7c42.inc> 

LIST 

:P17C42.INC Standard Header File, Version 1.03 Microchip Technology, Inc. 
LIST 



#define true 
#define false 
#define TRUE 
#define FALSE 

#define LSB 
#define MSB 



_INC equ 
_NO_INC egu 
_LOW equ 
_HIGH egu 



LIST 



1 



1 

#include <17c42 .mac> 



#include <17c42iir .mao 

LIST 

******************************************************************* 

PIC17C42 MACRO 

Macro For A Bi-Quad IIR Filter 

2nd order Direct Form (Transposed) Type 

Filter co-efficients B0 & B2 are assumed equal 

The difference equations for each cascade section is given by : 

Y(n) = B0*D(n) + Bl*D(n-l) + B2*D(n-2) 

D(n) = X(n) - Al*D(n-l) - A2*D(n-2) 
where X(n) = input sample, Y(n) = output of filter 
and Al, A2, B0, Bl, B2 are the Filter Co-efficients 

The above difference equations are only for 1 section of a 
2nd order Direct_Form II Filter structure (IIR) 

NOTE : 

It is possible to design the above structures 

such that the co-efficients B0 = B2 . If this is the 

case, 

Y(n) = B0*[D(n) + D(n-2)] = B2*[D(n) + D(n-2)] 
This way, one multiplication can be avoided 

If a 4th order filter is to be implemented, the output of 
the 1st structure should be input to the 2nd cascade section 
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00000001 



00029 
00030 
00031 
00032 
00033 
00034 
00035 
00036 
00037 
00038 
00039 
00040 
00041 
00042 
00043 
00044 
00045 
00046 
00047 
00048 
00049 
00050 
00051 
00052 
00053 
00054 
00055 
00056 
00057 
00058 
00059 
00060 
00061 
00062 
00063 
00064 
00065 
00066 
00067 
00068 
00069 
00070 
00071 
00072 
00073 
00074 
00075 
00076 
00077 
00078 
00079 
00080 
00081 
00082 
00083 
00084 
00085 
00086 
00087 
00088 
00089 
00090 
00091 
00092 
00093 
00094 



Timing (WORST CASE) s 

59+4*179 = 775 Cycles 
(194 uS a 16 Mhz) 
Program Memory : 
63 locations 

The sample filters are designed so that B0=B2 
This saves 1 multiplication 



B0_EQUALS_B2 equ 



TRUE 



**** + ****ir*************-k-k**ieir*****ir-k****** + ****-k**ir***+***** + + + **iri 
Parameters to BIQUAD Macro : 

Filter Constants Al, A2 , B0, Bl, B2 
& D(n), D(n-l), D(n-2), filter stage # 

BIQUAD MACRO Axl , Ax2 , BxO , Bxl , Dn , Dn_l , Dn_2 , stage 



; Compute Ax2*D(n-2) 
? 

MOVFP16 Dn_2,AARG ; D(n-2) = multiplier 

MOVFP16 Ax2 , BARG ; A2 = multiplicand 

call DblMult ; (ACCd,ACCc) = A2*D(n-2) 

Add product to output of 1st section 
Save result in 32 bit Accumulator 

ADD32 DPX,ACC 

Compute Al*D(n-l) 

MOVFP1 6 Dn_l , AARG ; AARG = D(n-2) = multiplier 

MOVFP16 Axl, BARG ; BARG = A2 = multiplicand 

call DblMult ; (ACCd,ACCc) = Al*D(n-l) 

Compute Al*D(n-l) + A2*D(n-2) + output of previous section 

multiplications already done, so simply perform a 32 bit add 
of previously obtained multiplication results 

ADD32 DPX,ACC ; ACC = Al*D(n-l) + A2*D(n-2) + (output of 1st section) 



; save the upper 16 bits of D(n) from the 32 bit accumulator 
; left shift the result by 1, to adjust the decimal point after 
; a Q15*Q15 multiplication 



rlcf 

rlcf 

movwf 

rlcf 

movwf 



ACC+Bl,w 
ACC+B2 , w 
Dn 

ACC+B3 , w 
Dn+Bl 



decimal adjust ( mult by 2) 



Compute B2 * [D(n) + D(n-2)] 
#if B0_EQUALS_B2 
ADD1 6 ACC Dn_2 , Dn , AARG 



; AARG = Dn + D(n-2) = multiplier 



M0VFP16 BxO, BARG ; BARG = A2 = multiplicand 

call DblMult ; (ACCd,ACCc) = B2* [D (n) +D (n-2 ) ] 

MOVPF32 DPX,ACC 

#else 

M0VFP16 BxO, BARG 
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00095 
00096 
00097 
00098 
00099 
00100 
00101 
00102 
00103 
00104 
00105 
00106 
00107 
00108 
00109 
00110 
00111 
00112 
00113 
00114 
00115 
00116 
00117 
00118 
00119 
00120 
00121 
00122 
00123 
00124 
00125 
00126 
00127 
00128 
00129 
00130 
00131 
00132 
00023 
00024 
00025 
00026 
00027 
00028 
00029 
00030 
00031 
00032 
00033 
00034 
00035 
00036 
00037 
00038 
00039 
00040 
00041 
00042 
00043 
00044 
00045 
00046 
00047 
00048 
00049 
00050 



M0VFP16 Dn , AARG 

call DblMult 
MOVPF32 DPX,ACC 

MOVFP16 Bx2 , BARG 
MOVFP16 Dn_2 , AARG 

call DblMult 
ADD32 DPX.ACC 



B0*D(n) 



t B2*D(n-2) 



#endif 

Shift down D(n-l) to D(n-2) after D(n-2) usage is no longer required. 
This way in the next iteration D(n-2) is equal to the present D(n-l) 

movf p Dn_l , AARG+B0 

movpf AARG+B0 , Dn_2 j Shift down D(n-l) 

movfp Dn_l+Bl , AARG+B1 

movpf AARG+B1 , Dn_2 +B1 ; AARG = D(n-l) = multiplier 



MOVFP16 Bxl , BARG 

call DblMult 



; BARG = Bl = multiplicand 

; (ACCd,ACCc) = Bl*D(n-l) 



( Compute Output Y = Bl*D(n-l) + B2*D(n-2) + B0*D(n) 

= Bl*D(n-l) + B0*[D(n) + D(n-2)] 
; Since all multiplications are already done, simply perform a 
; 32 bit addition 



ADD32 



DPX,ACC ; ACC = Bl*D(n-l) + B2*D(n-2) + B0*D (n) 



Shift down D(n) to D(n-l) so that in the next iteration, the new 
D(n-l) is the present D(n) 

MOV16 Dn,Dn_l ; Shift down D (n) to D(n-l) 

ENDM 

LIST 

******************************************************************* 

Second Order Direct Form IIR Filter 



In the code given below, a 4th order IIR Elliptic Lowpass Filter 
is implemented. Other order filters may be implemented by 
taking the following example code as a basis. 



Program: 
Revision Date: 



IIR. ASM 
1-13-97 



Compatibility with MPASMWIN 1.40 



******************************************************************* 

The specifications of the filter are : 
Sampling Frequency =2.0 Khz 

Filter Type = 4th Order Elliptic Lowpass Filter 



Lower Band Edge . 
Upper Band Edge 
Nominal Gain 
Nominal Ripple 0.01 



Bandl Band2 

600 Hz 
500 Hz 1 Khz 

1.0 0.0 
.05 



Maximum Ripple 0.00906 0.04601 
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00000000 



00000018 
0000001C 



00000020 
00000022 
00000024 

00000026 
00000028 



00051 
00052 
00053 
00054 
00055 
00056 
00057 
00058 
00059 
00060 
00061 
00062 
00063 
00064 
00065 
00066 
00067 
00068 
00069 
00070 
00071 
00072 
00073 
00074 
00075 
00076 
00077 
00078 
00079 
00080 
00081 
00082 
00083 
00084 
00085 
00086 
00087 
00088 
00089 
00090 
00091 
00092 
00093 
00094 
00095 
00096 
00097 
00098 
00099 
00100 
00101 
00102 
00103 
00104 
00105 
00106 
00107 
00108 
00109 
00110 
00111 
00112 
00113 
00114 
00115 
00116 



Ripple in dB 



0.07830 



-26.75 



The Filter Co-efficients for the above specifications 
of the filter are computed as follows : 



1st Section 



2nd Section 



All 
A12 
B10 
Bll 
B12 

A21 
A22 
B20 
B21 
B22 



-0.133331 

0.167145 

0.285431 

0.462921 

0.285431 

. 147827 
0.765900 
0.698273 
.499908 
0.698273 



Performance (WORST Case) : 



Cycles = #of Filter Stages*775 + 16 
= 2*775+16 = 1566 Cycles 
( 391 uSec) 
per each sample. Initialization 
time after reset is not counted 
Timing measured with B0_EQUALS_B2 
set to TRUE (see B I QUAD Macro for 
explanation) 



Program Memory 



= 16+ # of FilterStages * (BIQUAD Memory 
+ filter co-efficients) 
+ multiplier 

= 16+2* (63+5) +274 = 421 locations 
(excluding initialization) 

RAM usage = 48 file registers 

RAM usage/each additional stage = 16 file regs 



This time is less than 2 Khz (500 uSec) , 
which means real time filtering is possible 



CBLOCK 

BB0 , BB1 , BB2 , BB3 

ENDC 



CBLOCK 0x18 

DPX,DPX1,DPX2,DPX3 
AARG, AARG1 , BARG, BARG1 

ENDC 



; arithmetic accumulator 
; multiply arguments 



CBLOCK 



Dnl , Dnl_Hi 

Dnl_l, Dnl_l_Hi 

Dnl_2, Dnl_2_Hi 

Dn2, Dn2_Hi 

Dn2_l, Dn2_l_Hi 
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0000002A 


00117 


Dn2_ 


2 , Dn2_2_Hi 






00118 


ENDC 








00119 










00120 


CBLOCK 






0000002C 


00121 


All, 


All_Hi 




0000002E 


00122 


A12, 


A12_Hi 




00000030 


00123 


BIO, 


B10_Hi 




00000032 


00124 


Bll, 


Bll_Hi 


; 1st Section Filter Co-efficients 


00000034 


00125 


B12, 


B12_Hi 






00126 








00000036 


00127 


A21, 


A21_Hi 




00000038 


00128 


A22, 


A22_Hi 




0000003A 


00129 


B20, 


B20_Hi 




0000003C 


00130 


B21, 


B21_Hi 


; 2nd Section Filter Co-efficients 


0000003E 


00131 


B22, 


B22_Hi 






00132 


ENDC 








00133 










00134 


CBLOCK 






00000040 


00135 


X, XI 


; 16 bits of input stream 


00000042 


00136 


Y, Yl 


; 16 bits of filter output 




00137 








00000044 


00138 


ACC, 


ACC1, ACC2, ACC3 


,- 32 bit accumulator for computations 




00139 


ENDC 








00140 








00000002 


00141 FltStage .set 


2 




0000000A 


00142 NumCoeff equ 


(5*FltStage) 


; 5 Co-eff per stage 




00143 








00000001 


00144 


? equ 


1 






00145 








00000001 


00146 LPASS .set 


TRUE 




00000000 


00147 HPASS .set 


FALSE 




00000000 


00148 BPASS .set 


FALSE 


; select the desired filter type 


00000000 


00149 BSTOP .set 


FALSE 






00150 








00000001 


00151 SIGNED equ 


TRUE 


i Set This To 'TRUE' for signed 




00152 






; multiplication and 1 FALSE ' for unsigned 




00153 










00154 


*********++***+*+***+***+++****+****++*+********++****+**+********* 




00155 




Test Program For Low Pass Filter 




00156 


*+****+****+*******+*+**************+*+**+**++*+******+*****+*++**+ 




00157 








0000 


00158 


ORG 


0x0000 






00159 








0000 


00160 start 






0000 EOOD 


00161 


call 


initFilter 




0001 BOOO 


00162 


movlw 


0x00 




0002 0140 


00163 


movwf 


X 




0003 B07F 


00164 


movlw 


0x7 f 


; set initial Xn = X(0) = 0x7f00 


0004 0141 


00165 


movwf 


X+BB1 


; test for impulse response 




00166 








0005 


00167 NextPoint 






0005 E022 


00168 


call 


IIR_Filter 




0006 A442 


00169 


tlwt 


_LOW, Y 




0007 


00170 tracePoint 






0007 AE43 


00171 


tablwt 


_HIGH, 0, Y+BB1 




0008 0000 


00172 


nop 






0009 2940 


00173 


clrf 


X, F 


; set X(n) =0 , n <> 


OOOA 2941 


00174 


clrf 


X+BB1, F 


; for simulating an Impulse 


OOOB C005 


00175 


goto 


NextPoint 






00176 








OOOC COOC 


00177 self goto 


self 






00178 










00179 


******************************************************************* 




00180 , 








OOOD 


00181 initFilter 








00182 , 
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OOOD BOC1 
OOOE OlOD 
OOOF BOOl 
OOIO 010E 



00183 
00184 
00185 
00186 
00187 
00188 
00189 
00190 
00191 
00192 
00193 
00194 
00195 
00196 
00197 
00198 
00199 
00200 
00201 
00202 
00203 
00204 
00205 
00206 
00207 
00208 
00209 
00210 



At first read the Filter Co-efficients from Prog. Mem to Data RAM 
#if 



#endif 
#if 



#endif 
#if 



#endif 
#if 



#endif 



LPASS 
movlw 
movwf 
movlw 
movwf 

HPASS 
movlw 
movwf 
movlw 
movwf 

BPASS 
movlw 
movwf 
movlw 
movwf 

BSTOP 
movlw 
movwf 
movlw 
movwf 



LOW _coeff_lpass 
TBLPTRL 

HIGH _coeff_lpass 
TBLPTRH 



LOW _coeff_hpass 
TBLPTRL 

HIGH _coeff_hpass 
TBLPTRH 



LOW 

TBLPTRL 

HIGH 

TBLPTRH 



coef f_bpass 
coef f_bpass 



LOW _coeff_bstop 
TBLPTRL 

HIGH _coeff_bstop 
TBLPTRH 



0011 


B02C 


00211 


movlw 


All 


0012 


0101 


00212 


movwf 


FSR0 


0013 


8404 


00213 


bsf 


ALUSTA, FS0 


0014 


8D04 


00214 


bcf 


ALUSTA, FS1 j auto increment 






00215 










00216 


j Read Filter Co 


-efficients from Program Memory 






00217 






0015 


B0 OA 


00218 


movlw 


NumCoef f 


0016 


A92C 


00219 


tablrd 


_LOW,_INC,All ; garbage 


0017 




00220 


Next Coef f 




0017 


A000 


00221 


tlrd 


_LOW, INDF0 


0018 


AB00 


00222 


tablrd 


_HIGH,_INC, INDF0 


0019 


170A 


00223 


decf sz 


WREG, F 


001A 


C017 


00224 


goto 


NextCoef f 






00225 


t 








00226 


; Initialize "Dn"s to zero 






00227 






001B 


B020 


00228 


movlw 


Dnl 


001C 


0101 


00229 


movwf 


FSR0 


001D 


B00C 


00230 


movlw 


6*FltStage 


001E 




00231 


NextClr 




001E 


2900 


00232 


clrf 


INDF0, F 


001F 


170A 


00233 


decf sz 


WREG, F 


0020 


C01E 


00234 


goto 


NextClr 






00235 






0021 


0002 


00236 


return 








00237 










00238 


;**********************************************+***+**************** 






00239 




1st Cascade Section 






00240 








00241 






0022 




00242 


IIR Filter 





00243 
00244 
00245 
00246 
00247 
00248 



Compute D(n) = X(n) + Al*D(n-l) + A2*D(n-2) 

Since the filter constants are computated in Q15 format, 
X(n) must be multiplied by 2**15 and then added to the 
other terms . 



DS00540C-page 10 



© 1997 Microchip Technology Inc. 



AN540 



00249 
00250 



002C 7C24 
002D 7D25 



002E 7E2E 
002F 7P2F 



0030 EOAF 



0031 
0032 
0033 
0034 
0035 
0036 
0037 
0038 



6A18 
0F44 
6A19 
1145 
6A1A 
1146 
6A1B 
1147 



0039 7C22 
003A 7D23 



003B 7E2C 
003C 7F2D 



003D EOAF 



; Move Input to accumulator after proper scaling 



Scale the input X 



0022 


8804 


00251 


bcf 


ALUSTA,C 


0023 


1941 


00252 


rrcf 


X+BB1, F 


0024 


1940 


00253 


rrcf 


X, F 


0025 


290A 


00254 


clrf 


WREG, F 


0026 


190A 


00255 


rrcf 


WREG, F 


0027 


0145 


00256 


movwf 


ACC+BB1 


0028 


6A40 


00257 


movfp 


X,WREG 


0029 


0146 


00258 


movwf 


ACC+BB2 


002A 


6A41 


00259 


movfp 


X+BB1,WREG 


002B 


0147 


00260 


movwf 


ACC+BB3 



00261 
00262 
00263 
00264 
M 

M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
H 
H 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
H 
M 
M 
M 
M 
M 
M 
M 

M 
M 
M 
M 
M 
M 
M 
M 
M 



ACC = scaled input : X* (2**15) 
1st Biquad filter section 

BIQUAD Al 1 , A12 , B10 , Bll , Dnl , Dnl_l , Dnl_2 , 1 



Compute Ax2*D(n-2) 
MOVFP16 Dnl_2 , AARG 



; D(n-2) = multiplier 



MOVFP Dnl_2+B0,AARG+B0 ; move A (B0) to B(B0) 
MOVFP Dnl 2+Bl,AARG+Bl ; move A(B1) to B(B1) 



MOVFP1 6 A12,BARG 

MOVFP A12 +B0 , BARG+B0 
MOVFP A12 +B1 , BARG+B1 



call 



DblMult 



; A2 = multiplicand 

; move A (B0) to B(B0) 

; move A(B1) to B(B1) 

; (ACCd,ACCc) = A2*D(n-2) 



; Add product to output of 1st section 
,- Save result in 32 bit Accumulator 



ADD32 

MOVFP 

ADDWF 

MOVFP 

ADDWFC 

MOVFP 

ADDWFC 

MOVFP 

ADDWFC 



DPX,ACC 

DPX+B0,WREG 
ACC+B0, F 
DPX+B1 , WREG 
ACC+B1, F 
DPX+B2 , WREG 
ACC+B2, F 
DPX+B3 , WREG 
ACC+B3 , F 



Compute Al*D(n-l) 
MOVFP1 6 Dnl_l , AARG 



MOVFP 
MOVFP 



Dnl_l+B0 , AARG+B0 
Dnl 1+B1,AARG+B1 



MOVFP1 6 All.BARG 

MOVFP A11+B0,BARG+B0 
MOVFP A11+B1,BARG+B1 



get lowest byte of a into w 

add lowest byte of b, save in b(B0) 

get 2nd byte of a into w 

add 2nd byte of b, save in b(Bl) 

get 3rd byte of a into w 

add 3rd byte of b, save in b(B2) 

get 4th byte of a into w 

add 4th byte of b, save in b(B3) 



; AARG = D(n-2) = multiplier 

; move A(B0) to B(B0) 

; move A(B1) to B(B1) 

; BARG = A2 = multiplicand 



call 



DblMult 



A(B0) to B(B0) 
; move A(B1) to B(B1) 

; (ACCd,ACCc) = Al*D(n-l) 



Compute Al*D(n-l) + A2*D(n-2) + output of previous section 

multiplications already done, so simply perform a 32 bit add 
of previously obtained multiplication results 



ADD32 DPX,ACC 



; ACC = Al*D(n-l) + A2*D(n-2) + 
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M 






(output of 1st section) 


003E 


6A18 


M 


MOVFP 


DPX+B0,WREG 


get lowest byte of a into w 


003F 


0F44 


M 


ADDWF 


ACC+BO, F 


add lowest byte of b, save in b(B0) 


0040 


6A19 


M 


MOVFP 


DPX+B1 , WREG 


get 2nd byte of a into w 


0041 


1145 


M 


ADDWFC 


ACC+B1, F 


add 2nd byte of b, save in b(Bl) 


0042 


6A1A 


M 


MOVFP 


DPX+B2,WREG 


get 3rd byte of a into w 


0043 


1146 


M 


ADDWFC 


ACC+B2, F 


add 3rd byte of b, save in b(B2) 


0044 


6A1B 


M 


MOVFP 


DPX+B3,WREG 


get 4th byte of a into w 


004S 


1147 


M 


ADDWFC 


ACC+B3, F 


add 4th byte of b, save in b(B3) 



0046 1A45 

0047 1A46 

0048 0120 

0049 1A47 
004A 0121 



004B 6A24 
004C 0E20 
004D 011C 
004E 6A25 
004F 1021 
0050 011D 



0051 7E30 

0052 7F31 



0053 E0AF 



0054 5844 

0055 5945 

0056 5A46 

0057 5B47 



M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 



save the upper 16 bits of D(n) from the 32 bit accumulator 
left shift the result by 1, to adjust the decimal point after 
a Q15*Q15 multiplication 



rlcf ACC+B1 , w 

rlcf ACC+B2,w 

movwf Dnl 

rlcf ACC+B3,w 

movwf Dnl+Bl 



; decimal adjust ( mult by 2) 



Compute B2 * [D(n) + D(n-2)] 
#if B0_EQUALS_B2 
ADD16ACC Dnl 2 , Dnl , AARG 



AARG = Dn + D(n-2) = multiplier 



movf p 
addwf 
movwf 
movf p 
addwf c 
movwf 



Dnl_2+B0,WREG 
Dnl+BO, w 
AARG+B0 
Dnl_2+B1,WREG 
Dnl+Bl, w 
AARG+B1 



MOVFP16 B10 , BARG 



MOVFP 
MOVFP 



B10+B0, BARG+B0 
B10+B1.BARG+B1 



call DblMult 
MOVPF32 DPX,ACC 

MOVPF DPX+B0 , ACC+BO 

MOVPF DPX+B1 , ACC+B1 

MOVPF DPX+B2 , ACC+B2 

MOVPF DPX+B3 , ACC+B3 



#else 

MOVFP16 BIO, BARG 
MOVFP16 Dnl, AARG 
call DblMult 
MOVPF32 DPX,ACC 

MOVFP16 Bx2 , BARG 
MOVFP16 Dnl_2 , AARG 
call DblMult 
ADD32 DPX,ACC 



; BARG = A2 = multiplicand 

I move A (B0) to B(B0) 

; move A(B1) to B(B1) 

; (ACCd,ACCc) = B2* [D (n) +D (n-2 ) ] 



move A (B0) to B(B0) 

move A(B1) to B(B1) 

move A(B2) to B(B2) 

move A(B3) to B(B3) 



i B0*D(n) 



; B2*D(n-2) 



#endif 

Shift down D(n-l) to D(n-2) after D(n-2) usage is no longer required. 
This way in the next iteration D(n-2) is equal to the present D(n-l) 
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0058 7C22 

0059 5C24 
005A 7D23 
005B 5D25 



005C 7E32 
005D 7F33 



005E EOAF 



0067 6A20 

0068 0122 

0069 6A21 
006A 0123 



006B 7C2A 
006C 7D2B 



006D 7E38 
006E 7F39 



006F EOAF 



M 
M 
M 
M 
M 
M 
M 
M 
M 
H 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 



movfp Dnl_l, AARG+BO 

movpf AARG+BO , Dnl_2 

movfp Dnl_l+B1,AARG+Bl 

movpf AARG+Bl,Dnl_2+Bl 

MOVFP 16 Bll , BARG 



MOVFP 
MOVFP 



call 



B11+B0, BARG+B0 
Bll+Bl, BARG+B1 



DblMult 



; Shift down D(n-l) 

; AARG = D(n-l) = multiplier 

; BARG = Bl = multiplicand 

; move A (B0) to B(B0) 
; move A(B1) to B(B1) 

; (ACCd,ACCc) = Bl*D(n-l) 



Compute Output Y = Bl*D(n-l) + B2*D(n-2) + B0*D(n) 

= Bl*D(n-l) + B0*[D(n) + D(n-2)] 
Since all multiplications are already done, simply perform a 
32 bit addition 

ADD32 DPX,ACC ; ACC = Bl*D(n-l) + B2*D(n-2) + B0*D(n) 



005F 


6A18 


M 


MOVFP 


DPX+B0 , WREG 


get lowest byte of a into w 


0060 


0F44 


M 


ADDWF 


ACC+B0 , F 


add lowest byte of b, save in b(B0) 


0061 


6A19 


M 


MOVFP 


DPX+B1 , WR 


get 2nd byte of a into w 


0062 


1145 


M 


ADDWFC 


ACC+B1, F 


add 2nd byte of b, save in b(Bl) 


0063 


6A1A 


M 


MOVFP 


DPX+B2 , WREG 


get 3rd byte of a into w 


0064 


1146 


M 


ADDWFC 


ACC+B2 , F 


add 3rd byte of b, save in b(B2) 


0065 


6A1B 


M 


MOVFP 


DPX+B3 , WREG 


get 4th byte of a into w 


0066 


1147 


M 


ADDWFC 


ACC+B3 , F 


add 4th byte of b, save in b(B3) 



M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 

00265 
00266 
00267 
00268 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 



Shift down D(n) to D(n-l) so that in the next iteration, the new 
D(n-l) is the present D(n) 



MOV16 



Dnl,Dnl_l ; Shift down D(n) to D(n-l) 



MOVFP Dnl+BO , WREG 

MOVWF Dnl_l+B0 

MOVFP Dnl+Bl , WREG 

MOVWF Dnl 1+B1 



; get byte of a into w 
; move to b(B0) 
; get byte of a into w 
; move to b(Bl) 



2nd Biguad filter section 
BIQUAD A2 1 , A22 , B2 , B2 1 , Dn2 , Dn2_l , Dn2_2 , 2 



Compute Ax2*D(n-2) 

MOVFP 16 Dn2_2 , AARG 

MOVFP Dn2_2+B0, AARG+BO 
MOVFP Dn2_2+B1,AARG+Bl 

MOVFP 16 A22 , BARG 

MOVFP A2 2 +B0 , BARG+B0 
MOVFP A2 2 +B1 , BARG+B1 



call 



DblMult 



; D(n-2) = multiplier 

; move A (B0) to B(B0) 

; move A(B1) to B(B1) 

; A2 = multiplicand 

; move A (B0) to B(B0) 

; move A(B1) to B(B1) 

; (ACCd,ACCc) = A2*D(n-2) 



Add product to output of 1st section 
Save result in 32 bit Accumulator 



ADD32 



DPX,ACC 
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0070 
0071 
0072 
0073 
0074 
0075 
0076 
0077 



6A18 
0F44 
6A19 
1145 
6A1A 
1146 
6A1B 
1147 



0078 7C28 

0079 7D29 



007A 7E36 
007B 7F37 



007C EOAF 



007D 
007E 
007F 
0080 
0081 
0082 
0083 
0084 



6A18 
0F44 
6A19 
1145 
6A1A 
1146 
6A1B 
1147 



0085 1A45 

0086 1A46 

0087 0126 

0088 1A47 

0089 0127 



008A 6A2A 
008B 0E26 
008C 011C 
008D 6A2B 
008E 1027 
008F 011D 



M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 



MOVFP 

ADDWF 

MOVFP 

ADDWFC 

MOVFP 

ADDWFC 

MOVFP 

ADDWFC 



DPX+BO , WREG 
ACC+BO , F 
DPX+B1 , WREG 
ACC+B1, F 
DPX+B2 , WREG 
ACC+B2 , F 
DPX+B3 , WREG 
ACC+B3 , F 



Compute Al*D(n-l) 



MOVFP16 Dn2 1,AARG 



get lowest byte of a into w 

add lowest byte of b, save in b(B0) 

get 2nd byte of a into w 

add 2nd byte of b, save in b(Bl) 

get 3rd byte of a into w 

add 3rd byte of b, save in b(B2) 

get 4th byte of a into w 

add 4th byte of b, save in b(B3) 



AARG = D(n-2) = multiplier 



MOVFP 
MOVFP 



Dn2_l+B0,AARG+B0 
Dn2_l+Bl , AARG+B1 



MOVFP16 A2 1 , BARG 



MOVFP 
MOVFP 

call 



A21+B0, BARG+B0 
A21+B1.BARG+B1 

DblMult 



; move A (B0) to B(B0) 

; move A(B1) to B(B1) 

; BARG = A2 = multiplicand 

; move A (B0) to B(B0) 

; move A(B1) to B(B1) 

; (ACCd,ACCc) = Al*D(n-l) 



Compute Al*D(n-l) + A2*D(n-2) + output of previous section 

multiplications already done, so simply perform a 32 bit add 
of previously obtained multiplication results 



ACC = Al*D(n-l) + A2*D(n-2) + 

(output of 1st section) 
get lowest byte of a into w 
add lowest byte of b, save in b(B0) 
get 2nd byte of a into w 



ADD32 


DPX , ACC 


MOVFP 


DPX+BO, WREG 


ADDWF 


ACC+BO, F 


MOVFP 


DPX+B1 , WREG 


ADDWFC 


ACC+B1 , F 


MOVFP 


DPX+B2 , WREG 


ADDWFC 


ACC+B2 , F 


MOVFP 


DPX+B3 , WREG 


ADDWFC 


ACC+B3 , F 



add 2nd byte of 
get 3rd byte of 
add 3rd byte of 
get 4th byte of 



save in b (Bl) 
into w 

save in b(B2) 
into w 



add 4th byte of b, save in b(B3) 



save the upper 16 bits of D (n) from the 32 bit accumulator 
left shift the result by 1, to adjust the decimal point after 
a Q15*Q15 multiplication 



rlcf ACC+Bl,w 

rlcf ACC+B2 , w 

movwf Dn2 

rlcf ACC+B3 , w 

movwf Dn2 +B1 



decimal adjust ( mult by 2) 



Compute B2 * [D(n) + D(n-2)] 

#if B0_EQUALS_B2 

ADD 16 ACC Dn2_2 , Dn2 , AARG 

movf p Dn2_2 +B0 , WREG 

addwf Dn2+B0,w 

movwf AARG+B0 

movf p Dn2_2 +B1 , WREG 

addwf c Dn2+Bl,w 

movwf AARG+B1 

MOVFP16 B2 , BARG ; 



AARG = Dn + D(n-2) = multiplier 



BARG = A2 = multiplicand 
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0090 7E3A 

0091 7F3B 



0092 EOAF 



0093 5844 

0094 5945 

0095 5A46 

0096 5B47 



0097 7C28 

0098 5C2A 

0099 7D29 
009A 5D2B 



009B 7E3C 
009C 7F3D 



009D EOAF 



009E 
009F 
00 AO 
00A1 
00A2 
00 A3 
00A4 
00A5 



6A18 
0F44 
6A19 
1145 
6A1A 
1146 
6A1B 
1147 



MOVFP 
MOVFP 



B20+B0,BARG+B0 
B20+B1, BARG+B1 



call DblMult 
MOVPF32 DPX,ACC 

MOVPF DPX+BO , ACC+BO 

MOVPF DPX+B1,ACC+B1 

MOVPF DPX+B2 , ACC+B2 

MOVPF DPX+B3 , ACC+B3 



#else 

MOVFP16 B2 , BARG 
MOVFP16 Dn2 , AARG 
call DblMult 
MOVPF3 2 DPX,ACC 

MOVFP16 Bx2,BARG 
MOVFP16 Dn2_2 , AARG 
call DblMult 
ADD32 DPX,ACC 

#endif 



move A(B0) to B(B0) 
move A(B1) to B(B1) 

(ACCd.ACCc) = B2* [D(n)+D(n-2) ] 



move A(B0) to B(B0) 

move A(B1) to B(B1) 

move A(B2) to B(B2) 

move A(B3) to B(B3) 



B0*D(n) 



B2*D(n-2) 



Shift down D(n-l) to D(n-2) after D(n-2) usage is no longer required. 
This way in the next iteration D(n-2) is equal to the present D(n-l) 



movfp Dn2_l , AARG+B0 

movpf AARG+B0 , Dn2_2 

movfp Dn2_l+Bl , AARG+B1 

movpf AARG+B1 , Dn2_2+Bl 

MOVFP16 B2 1 , BARG 



MOVFP 
MOVFP 



call 



B21+B0,BARG+B0 
B21+B1,BARG+B1 



DblMult 



; Shift down D(n-l) 

; AARG = D(n-l) = multiplier 

; BARG = Bl = multiplicand 

; move A(B0) to B(B0) 
; move A(B1) to B(B1) 

; (ACCd,ACCc) = Bl*D(n-l) 



Compute Output Y = Bl*D(n-l) + B2*D(n-2) + B0*D(n) 

= Bl*D(n-l) + B0*[D(n) + D(n-2)] 
Since all multiplications are already done, simply perform a 
32 bit addition 



ADD32 DPX,ACC 



MOVFP 

ADDWF 

MOVFP 

ADDWFC 

MOVFP 

ADDWFC 

MOVFP 

ADDWFC 



DPX+BO, WREG 
ACC+BO, F 
DPX+B1,WREG 
ACC+B1 , F 
DPX+B2 , WREG 
ACC+B2 , F 
DPX+B3 , WREG 
ACC+B3, F 



ACC = Bl*D(n-l) + B2*D(n-2) + B0*D(n) 

get lowest byte of a into w 

add lowest byte of b, save in b(B0) 

get 2nd byte of a into w 

add 2nd byte of b, save in b(Bl) 

get 3rd byte of a into w 

add 3rd byte of b, save in b(B2) 

get 4th byte of a into w 

add 4th byte of b, save in b(B3) 



Shift down D(n) to D(n-l) so that in the next iteration, the new 
D(n-l) is the present D(n) 



00A6 6A26 
00A7 0128 



M 
M 
M 
M 
M 
M 

M MOVFP 
M MOVWF 



MOV16 Dn2,Dn2_l 



Dn2 +B0 , WREG 
Dn2 1+B0 



; Shift down D(n) to D(n-l) 

; get byte of a into w 
; move to b(B0) 
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00A8 6A27 
00A9 0129 



00AA 6A46 
OOAB 0142 
00 AC 6A47 
OOAD 0143 



OOAE 0002 



01C1 

01C1 1111 
01C2 EA9B 
01C3 2489 
01C4 3B41 
01C5 2489 

01C6 ED14 
01C7 9DF7 
01C8 5961 
01C9 3FFD 
01CA 5961 



M 
M 
M 

M 

00269 
00270 
00271 
00272 
00273 
00274 
00275 
M 
M 
M 
M 
M 
M 

00276 
00277 
00278 
00279 
00280 
00281 
00282 
00178 
00283 
00284 
00285 
00286 
00287 
00288 
00289 
00290 
00291 
00292 
00293 
00294 
00295 
00296 
00297 
00298 
00299 
00300 
00301 
00302 
00303 
00304 
00305 
00306 
00307 
00308 
00309 
00310 
00311 
00312 
00313 
00314 
00315 
00316 
00317 
00318 
00319 
00320 
00321 
00322 
00323 



MOVFP Dn2 +B1 , WREG 
MOVWF Dn2 1+B1 



j get byte of a into w 
; move to b(Bl) 



The filter output is now computed 

Save the Upper 16 Bits of 3 2 bit Accumulator into Y after 
adjusting the decimal point 



MOV16 



ACC+BB2 , Y 



MOVFP ACC+BB2+B0,WREG 

MOVWF Y+B0 

MOVFP ACC+BB2+B1,WREG 

MOVWF Y+Bl 



return 



; get byte of a into w 
; move to b(B0) 
; get byte of a into w 
; move to b(Bl) 



; Output Y(n) computed 



. ******************************************************************* 

; Set SIGNED /UNSIGNED Flag Before Including 17c42MPY.mac 

#include <17c42MPY .mao 

LIST 

******************************************************************* 

Low Pass Filter Co-efficients 



ELLIPTIC 



LOWPASS FILTER 



FILTER ORDER = 4 
Sampling frequency = 

BAND 1 



LOWER BAND EDGE 
UPPER BAND EDGE 
NOMINAL GAIN 
NOMINAL RIPPLE 
MAXIMUM RIPPLE 
RIPPLE IN dB 

I A(I,1) 

1 -.133331 
.147827 



.00000 
.50000 
1.00000 
.01000 
.00906 
.07830 

MI, 2) 

.167145 
.765900 



2.000 KiloHertz 
BAND 2 



.60000 
1.00000 
.00000 
.05000 
.04601 
-26.74235 



B(I,0) 

.285431 
.698273 



B(I,1) 

.462921 
.499908 



B(I,2) 

.285431 
.698273 



_coeff_lpass ; co-efficients for 1st Cascade section 

data 4369 ; -All 

data -5477 ; -A12 

data 9353 ; B10 

data 15169 ; Bll 

data 9353 ; B12 

; co-efficients for 2nd Cascade section 

data -4844 ; -A21 

data -25097 ; -A22 

data 22881 ; B20 

data 16381 ; B21 

data 22881 ; B22 

******************************************************************* 
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00324 






00325 






00326 






00327 






00328 






00329 






00330 






00331 






00332 






00333 






00334 






00335 






00336 






00337 






00338 






00339 






00340 






00341 






00342 






00343 






00344 






00345 






00346 






00347 






00348 






00349 


01CB 




00350 


OICB 


DC8F 


00351 


OICC 


E6F5 


00352 


OICD 


2079 


00353 


OICE 


CB57 


003 54 


OICF 


2079 


003 55 






003 56 


OlDO 


0C12 


00357 


01D1 


9C1C 


00358 


01D2 


56DE 


00359 


01D3 


C1D0 


00360 


01D4 


56DE 


00361 






00362 






00363 






00364 






00365 






00366 






00367 






00368 






00369 






00370 






00371 






00372 






00373 






00374 






00375 






00376 






00377 






00378 






00379 






00380 






00381 






003 82 






00383 






00384 






00385 






00386 






00387 






00388 






00389 



High Pass Filter Co-efficients 

ELLIPTIC HIGHPASS FILTER 
FILTER ORDER = 4 

Sampling frequency = 2.000 KiloHertz 



LOWER BAND EDGE 
UPPER BAND EDGE 
NOMINAL GAIN 
NOMINAL RIPPLE 
MAXIMUM RIPPLE 
RIPPLE IN dB 



A(I,1) 

.276886 
■ .094299 



BAND 1 



.00000 
.40000 
.00000 
.04000 
.03368 
-29 .45335 

A(I,2) 

.195648 
.780396 



BAND 2 



.50000 
1.00000 
1.00000 
.02000 
.01686 
.14526 



B(I,0) 

.253677 
.678650 



B(I,1) 

• .411407 

• .485840 



B(I,2) 

.253677 
.678650 



_coef f_hpass 
data 
data 
data 
data 
data 

data 
data 
data 
data 
data 



; co-efficients for 1st Cascade section 

-9073 ; -All 

-6411 ; -A12 

8313 ; BIO 

-13481 ; Bll 

8313 ; B12 

; co-efficients for 2nd Cascade section 

3090 ; -A21 

-25572 ; -A22 

22238 ; B20 

-15920 ; B21 

22238 ; B22 



it****************************************************************** 

*******************************************^ 
Band Pass Filter Co-efficients 



ELLIPTIC 



BANDPASS FILTER 



FILTER ORDER = 4 
Sampling frequency = 

BAND 1 



2.000 KiloHertz 

BAND 2 BAND 3 



LOWER BAND EDGE 
UPPER BAND EDGE 
NOMINAL GAIN 
NOMINAL RIPPLE 
MAXIMUM RIPPLE 
RIPPLE IN dB 



A(I,1) 



- . 936676 
. 936707 



.00000 
.10000 
.00000 
.05000 
.03644 
-28.76779 



A(I,2) 

. 550568 
. 550568 



.30000 
.70000 
1.00000 
.05000 
.03867 
.32956 



B(I,0) 

.444000 
.615540 



. 90000 
1.00000 
.00000 
. 05000 
. 03641 
-28.77647 



B(I,1) 



- .865173 
1 .199402 



B(I,2) 

.444000 
.615540 
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00390 














01D5 




00391 


_coef f_ 


bpass 


) CO 


■efficients for 1st Cascade section 


01D5 


3BF2 


00392 




data 


30693/2 






-All 


01D6 


DCC4 


00393 




data 


-18041/2 






-A12 


01D7 


1C6A 


00394 




data 


14549/2 






B10 


01D8 


C8A1 


00395 




data 


-28350/2 






Bll 


01D9 


1C6A 


00396 




data 


14549/2 






B12 






00397 






; co-efficients for 2nd Cascade section 


01DA 


C40D 


00398 




data 


-30694/2 






-A21 


OlDB 


DCC4 


00399 




data 


-18041/2 






-A22 


OlDC 


2765 


00400 




data 


20170/2 






B20 


OIDD 


4CC3 


00401 




data 


39302/2 






B21 


OlDE 


2765 


00402 
00403 




data 


20170/2 






B22 






00404 


****** 


************************************************************* 






00405 


















00406 


****** 


************************************************************* 






00407 


Banc 


Stop Filter Co-efficients 










00408 


















00409 


















00410 






ELLIPTIC 




BANDSTOP FILTER 






00411 


















00412 






FILTER ORDER = 


4 








00413 






Sampling frequency 




2.000 KiloHertz 






00414 


















00415 






BAND 1 






BAND 2 BAND 3 






00416 


















00417 


















00418 


LOWER 


BAND EDGE .00000 






.45000 .70000 






00419 


UPPER 


BAND EDGE 


.30000 






.55000 1.00000 






00420 


NOMINAL GAIN 


1.00000 






.00000 1.00000 






00421 


NOMINAL RIPPLE 


.05000 






.05000 .05000 






00422 


MAXIMUM RIPPLE 


.03516 






.03241 .03517 






00423 


RIPPLE 


IN dB 


.30015 






-29.78523 .30027 






00424 


















00425 


















00426 


I 


A(1, 1) 


A(I,2) 






B(I,0) B(I,1) B(I,2) 






00427 


















00428 


1 


.74942C 


. 583282 






.392685 .087936 .392685 






00429 


2 


-.749390 .583282 






1.210022 -.270935 1.210022 






00430 














OlDF 




00431 


coef f_bstop 


; co- 


efficients for 1st Cascade section 


OlDF 


DO OA 


00432 




data 


-24557/2 




All 


OlEO 


DAAC 


00433 




data 


-19113/2 




A12 


01E1 


1922 


00434 




data 


12868/2 


B10 


01E2 


05 AO 


00435 




data 


2881/2 


Bll 


01E3 


1922 


00436 




data 


12868/2 


B12 






00437 






; co-efficients for 2nd Cascade section 


01E4 


2FF6 


00438 




data 


24556/2 




A21 


01E5 


DAAC 


00439 




data 


-19113/2 




A22 


01E6 


4D71 


00440 




data 


39650/2 


B20 


01E7 


EEA9 


00441 




data 


-8878/2 


B21 


01E8 


4D71 


00442 




data 


39650/2 


B22 






00443 


















00444 


******************************************************************* 






00445 


















00446 




END 
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MEMORY USAGE MAP ('X' = Used, = Unused) 

0000 : XXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXX 

0040 : XXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXX 

0080 : XXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXX 

00C0 : XXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXX 

0100 : XXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXX 

0140 : XXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXX 

0180 : XXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXX 

01C0 : XXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXX xxxxxxxxx 

All other memory blocks unused. 

Program Memory Words Used: 489 



Errors : 

Warnings : reported, suppressed 

Messages : reported, suppressed 
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Note the following details of the code protection feature on PICmicro® MCUs. 

The PICmicro family meets the specifications contained in the Microchip Data Sheet. 

Microchip believes that its family of PICmicro microcontrollers is one of the most secure products of its kind on the market today, 
when used in the intended manner and under normal conditions. 

There are dishonest and possibly illegal methods used to breach the code protection feature. All of these methods, to our knowl- 
edge, require using the PICmicro microcontroller in a manner outside the operating specifications contained in the data sheet. 
The person doing so may be engaged in theft of intellectual property. 

Microchip is willing to work with the customer who is concerned about the integrity of their code. 

Neither Microchip nor any other semiconductor manufacturer can guarantee the security of their code. Code protection does not 
mean that we are guaranteeing the product as "unbreakable". 

Code protection is constantly evolving. We at Microchip are committed to continuously improving the code protection features of 
our product. 

If you have any further questions about this matter, please contact the local sales office nearest to you. 



Information contained in this publication regarding device 
applications and the like is intended through suggestion only 
and may be superseded by updates. It is your responsibility to 
ensure that your application meets with your specifications. 
No representation or warranty is given and no liability is 
assumed by Microchip Technology Incorporated with respect 
to the accuracy or use of such information, or infringement of 
patents or other intellectual property rights arising from such 
use or otherwise. Use of Microchip's products as critical com- 
ponents in life support systems is not authorized except with 
express written approval by Microchip. No licenses are con- 
veyed, implicitly or otherwise, under any intellectual property 
rights. 



Trademarks 

The Microchip name and logo, the Microchip logo, FilterLab, 
KeeLoq, microlD, MPLAB, PIC, PICmicro, PICMASTER, 
PICSTART, PRO MATE, SEEVAL and The Embedded Control 
Solutions Company are registered trademarks of Microchip Tech- 
nology Incorporated in the U.S.A. and other countries. 

dsPIC, ECONOMONITOR, FanSense, Flex ROM, fuzzyLAB, 
In-Circuit Serial Programming, ICSP, ICEPIC, microPort, 
Migratable Memory, MPASM, MPLIB, MPLINK, MPSIM, 
MXDEV, PICC, PICDEM, PICDEM.net, rfPIC, Select Mode 
and Total Endurance are trademarks of Microchip Technology 
Incorporated in the U.S.A. 

Serialized Quick Turn Programming (SQTP) is a service mark 
of Microchip Technology Incorporated in the U.S.A. 

All other trademarks mentioned herein are property of their 
respective companies. 

© 2002, Microchip Technology Incorporated, Printed in the 
U.S.A., All Rights Reserved. 
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ISO 9001 / QS-9000 
REGISTERED FIRM 



Microchip received QS-9000 quality system 
certification for its worldwide headquarters, 
design and wafer fabrication facilities in 
Chandler and Tempe, Arizona in July 1999. The 
Company's quality system processes and 
procedures are QS-9000 compliant for its 
PICmicro" 8-bit MCUs, KeeLoq" code hopping 
devices, Serial EEPROMs and microperipheral 
products. In addition, Microchip's quality 
system for the design and manufacture of 
development systems is ISO 9001 certified. 
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